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ABSTRACT 

We present our investigation of pre-main-sequence (PMS) stellar populations in the Large Magellanic Cloud 
(LMC) from imaging with Hubble Space Telescope WFPC2 camera. Our targets of interest are four star- 
forming regions located at the periphery of the super-giant shell LMC 4 (Shapley Constellation III). The 
PMS stellar content of the regions is revealed through the differential Hess diagrams and the observed color- 
magnitude diagrams (CMDs). Further statistical analysis of stellar distributions along cross-sections of the 
faint part of the CMDs allowed the quantitative assessment of the PMS stars census, and the isolation of faint 
PMS stars as the true low-mass stellar members of the regions. These distributions are found to be well rep- 
resented by a double Gaussian function, the first component of which represents the main-sequence field stars 
and the second the native PMS stars of each region. Based on this result, a cluster membership probability 
was assigned to each PMS star according to its CMD position. The higher extinction in the region LH 88 did 
not allow the unambiguous identification of its native stellar population. The CMD distributions of the PMS 
stars with the highest membership probability in the regions LH 60, LH 63 and LH 72 exhibit an extraordinary 
similarity among the regions, suggesting that these stars share common characteristics, as well as common 
recent star formation history. Considering that the regions are located at different areas of the edge of LMC 4, 
this finding suggests that star formation along the super-giant shell may have occurred almost simultaneously. 

Subject headings: Magellanic Clouds - Hll regions - Hertzsprung-Russell and C-M diagrams - open clusters 
and associations: individual (LH 60, LH 63, LH 72, LH 88) - stars: formation - stars: 
pre-main-sequence 



1. INTRODUCTION 

Considerable star formation in the Large Magellanic Cloud 
(LMC) takes place in interstellar shells, ranging from small 
bubbles to large super-g iant shells (see , e.g., |Chu 2009). 
Among the latter, LMC 4 (Meabum 1980) is the largest super- 
giant shell in the Local Group, encompassing a remarkable 
Hi ca vity of diameter ^ L9 kpc dDopita. Mat hewson. & Ford 
119851) . Centered on Shapley Constellation III (Shapley. 1951.) 
the area of LMC 4 compris es over 500 clu sters, associa- 
tions and emission nebulae (IBica et alJll999h . In the Milky 
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Way, the youngest stellar associations with r < 10 Myr are 
surrounded by bright Hll regions, and comprise high-mass 
main sequence (MS) stars as well as intermediate- and low- 
mass pre-main sequence (PMS) stars. While the high-mass 
(OB-type) MS and intemediate-mass (Herbig Ae/Be) PMS 
populations are the direct signature of the youthfulness of 
their hosting associations, the low-mass PMS stars preserve 
a record of the complete recent star formation history of the 
region over long periods, since their evolution is extremely 
slow and can last up to many tens of MyiQ Bearing this in 
mind, we undertake a research project that aims at a compre- 
hensive study of the stellar populations in LMC star-forming 
regions, with emphasis on the recent star formation history 
in the vicinity of LMC 4 as recorded in the low-mass PMS 
populations. 

The existence of such stars in star-forming regions of the 
LMC became recently known thanks to the angular reso- 
lution and wide-field coverage provided by Hubble Space 
Telescope (HST). Archival images taken with the Wide-Field 
Planetary Came ra 2 (WF PC2) of the young LMC associa- 
tion LH 52 (Lu cke & Hod ge 1970), located at the northeast- 
em edge of LMC 4, revealed for the first time that low-mass 
PMS stars can be directly identified in the color-magnitude 
diagram (CMD) from photometry in V- and /-equivalent 
filters (Gouliermis, Brandner, & Henning 2006). While the 
WFPC2 images of LH 52 provided the first proof of the exis- 
tence of such young stars in LMC star-forming regions, they 
were not deep enough to allow a statistically sound investi- 
gation of these stars. Subsequent deep photometry with the 

Typical contraction time for a 1 Mq star is 50 Myr and for a 0.5 M0 
star 200 Myr ( Karttunen et al. 2 00% . 
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TABLE 1 

Description of the WFPC2 observations within our HST Program GO- 1 1547. 



Association Visit 


P A 
K./\. 


ueci. 


Data set 


r liter 




Lxposure 


Observation 


l^r^icia^ i>o 






lilename 






time (s) 


Date 


LH63 1 


05''27"M7!91 


-67°25'43'.'0 


ubOhOlO 


F300W 




1 X 10, 6 X 160 


2008 Jul 15 










F450W 




1 X 10, 6 X 160 


2008 Jul 15 










F555W 




4 X 300, 2 X 350 


2008 Jul 15 










F656N 


1 


X 300, 1 X 260, 1 X 800 


2008 Jul 15 










F814W 




4 X 260, 2 X 300 


2008 Jul 15 


2 


05'^28™05!60 


-67°25'35'.'l 


ub0h020 


F555W 




4 X 300, 2 X 350 


2008 Jul 19 










F814W 




4 X 260, 2 X 300 


2008 Jul 19 


LH 60 3 


051^ 27™ IK 89 


-67°27'27'.'3 


ubOh030 


F300W 




1 X 10, 6 X 160 


2008 Jul 18 










F450W 




1 X 10, 6 X 160 


2008 Jul 18 










F555W 




4 X 300, 2 X 350 


2008 Jul 18 










F656N 


1 


X 300, 1 X 260, 1 X 800 


2008 Jul 18 










F814W 




4 X 260, 2 X 300 


2008 Jul 18 


4 


05'^27'"28?56 


~67°27'19'.'4 


ub0h040 


F555W 




4 X 300, 2 X 350 


2008 Jul 19 










F814W 




4 X 260, 2 X 300 


2008 Jul 19 


LH72 5 


05'^32™26!03 


-66°28'15'.'8 


ubOhOSO 


F300W 




1 X 10, 6 X 160 


2008 Jul 21 










F450W 




1 X 10, 6 X 160 


2008 Jul 21 










F555W 




4 X 300, 2 X 350 


2008 Jul 21 










F656N 


1 


X 300, 1 X 260, 1 X 800 


2008 Jul 21 










F814W 




4 X 260, 2 X 300 


2008 Jul 21 


6 


05''32™08!59 


-66°28'24'.'0 


ubOhOeO 


F555W 




4 X 300, 2 X 350 


2008 Jul 14 










F814W 




4 X 260, 2 X 300 


2008 Jul 14 


LH 88 7 


05'^35™55!13 


-67°34'39'.'9 


ub0h570 


F300W 




1 X 10, 6 X 160 


2008 Sep 1 










F450W 




1 X 10, 6 X 160 


2008 Sep 1 










F555W 




4 X 300, 2 X 350 


2008 Sep 1 










F656N 


1 


X 300, 1 X 260, 1 X 800 


2008 Sep 1 










F814W 




4 X 260, 2 X 300 


2008 Sep 1 


8 


05'^35™37?72 


-67°34'50'.'5 


ubOhOSO 


F555W 




4 X 300, 2 X 350 


2008 Jul 9 










F814W 




4 X 260, 2 X 300 


2008 Jul 9 


Field 1 


05'^41™38!52 


-68°17'01'.'8 


u9e5550 


F555W 




2 X 500 


2008 Feb 21 










F814W 




2 X 500 


2008 Feb 21 


Field 2 


05h24™04!72 


-67°21'57'.'7 


u9px070 


F555W 




2 X 500 


2007 May 20 










F814W 




2 X 500 


2007 May 20 



Note. — Observations of Fields 1 and 2, for the assessment of the general LMC field population, are observed within HST 
programs GO-10583 (PI: C. Stubbs) and GO-10903 (PI: A. Rest) respectively. These data were retrieved from the HST Data 
Archive. 



Advanced Camera for Surveys (ACS) in the same filters of 
the association LH 95, located at the northwest periphery of 
LMC 4, revealed an outstanding sample of more than 2,500 
low-mass PMS stars down to the (model-dependent) mass- 
limit of ^ 0.2 Mp,, the smallest stellar mass ever observed in 
another galaxy (iGouliermis et al.|[2007h . 

These data are characterized by an unprecedented com- 
pleteness in their photometry due to their deepness. They 
provided us, thus, a unique sample of PMS stars, which we 
utilized to address in detail the formation of a young stel- 
lar cluster in another galaxy. Specifically, we addressed the 
initial mass function (IMF) of LH 95 down t o the sub-solar 
regime jDa Rio. Gouliermis. & Henning|20()9l) and the age of 
the cluster, as well as the duration of the star formation pro- 
cess in it, via the detection of an a ge-spread among its low- 
mass PMS stars (Da Rio. Goulierm is. & Genn aro 2010). The 
aforementioned investigations provide a unique insight of the 
faintest stars in two star-forming regions within LMC 4. How- 
ever, a more complete analysis of PMS populations in this 
area requires the collection of data on a larger sample of such 
objects. Within our HST program GO-11547 we obtained 
multi-band imaging of four additional star-forming regions 
located at the periphery of LMC 4 in order to accomplish a 
thorough characterization of PMS, MS and evolved stars in 
their vicinity, study their formation process, determine the 
stellar IMF and its variations in different environments, and 
provide an ample set of physical parameters of all observed 
stellar types in these regions. 



In this first paper we present the photometry, identify the 
low-mass PMS stellar populations in the observed regions, 
and determine their cluster membership via their CMD po- 
sitions. Specifically, in §|2]we present the observations of this 
program and describe the photometric process, as well as its 
accuracy and completeness. Our measurements of the inter- 
stellar reddening towards the targets, based on the photometry 
of the brighter MS stars, are presented in § [3] We present the 
observed CMDs in §|4l where we also discuss the CMD of the 
general LMC field. The identification and a qualitative as- 
sessment of the PMS stars in the observed regions takes place 
in §|5] where the Hess diagrams and a Monte Carlo technique 
for the statistical subtraction of field stars from the observed 
CMDs are described. In §|6]we present the quantitative deter- 
mination of the true PMS populations of the observed regions 
from the stellar distributions along cross-sections of the ob- 
served CMDs. We discuss our findings in § [T] Concluding 
remarks and our plans for subsequent analyses are given in 

m 

2. OBSERVATIONS AND PHOTOMETRY 

2.1. Observations 

Our investigation is based on the analysis of data collected 
within our HST program GO-11547 (PI: D. GouHermis) that 
imaged young stellar clusters located at the apparent edge of 
the super-giant shell LMC 4. Our targets are the young asso- 
ciations LH 60, LH 63, LH 72 and LH 88 ( Lucke & Hodgg 
[1970.) . embedded in the bright emission nebulae DEM L 201, 
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Fig. 1 . — Maps of the regions observed. The general area of the super-giant shell LMC 4 is shown in the center from a scanned UK Schmidt telescope plate 
observed in the R band, available from the DSS. The extracted subregions, shown on the left and right of the figure, cover about 10' X 7' around the targets. 
The observed WFPC2 pointings are overlaid and annotated according to the con'esponding visit number (see also Table [T]. The scale indicator corresponds to 
the fields-of-view of the individual regions around the observed clusters. The field-of-view of the image of the general LMC 4 area is ~ 1.7° X 2.2° (about 
1.5 kpc X 2 kpc). 



228 and 241 dPavies. Elliott. & Meaburnlll976l). coinciding 
with the Hll r egions LHA 120-N 51A,C, N 55A and N 59C 
(lHenizeiri956l) . respectively. Maps of the general region of 
LMC 4 and the immediate areas of the targets are shown in 
Fig. [U As such star-forming regions extend to sizes up to 
a few 10 pc (1' ~ 15 pc at the distance of the LMC), our 
investigation requires the best combination of high angular 
resolution, to overcome crowding effects for faint stars, and 
wide field-of-view to effectively sample each region. During 
Cycle 16 the most appropriate instrument for this task was 
WFPC2. To achieve large spatial coverage each cluster was 
observed with two overlapping pointings. The first pointing 
was centered on the brightest part of the associated Hll region 
and WFPC2 imaging was performed in the broad-band filters 
F300W, F450W, F555W, F814W and the narrow-band filter 
F656N, roughly equivalent to standard U, B, V, I and Ha, 
respectively. We refer to each of these first group of multi- 
band pointings on the Hll regions as Pointing 1. The second 
WFPC2 pointings of each cluster, referred to as Pointing 2, 
are centered off the bright nebula and were observed only in 
the F555W and F814W filters. 

Since we wish to eliminate confusion effects, especially for 
the detection of the low-mass PMS populations, all observa- 
tions are performed with a suitable dithering, using a sub- 
pixel pointing pattern, in order to improve the spatial resolu- 
tion. Furthermore, long exposures were split into intervals to 
aid in the removal of cosmic rays without any significant cost 
in the S/N of the final combined image. Short exposures in 
the F300W and F450W filters are used to obtain photometry 
for the brightest and youngest stars in these regions. Previous 
ground-based photometry is also available for these stars (e.g . , 
iHill. Madore. & Freedman|[T994t lOlsen. Kim. & Bussll2001h . 
The consideration of the guide star acquisition, reacquisition 
of the telescope, and the instrument overhead times (e.g., 
change of filter, CCD readout, dithering) resulted in one set of 
four orbits for each Pointing 1 and two orbits for each Point- 



ing 2. The orbits for each pointing were grouped as a visit, 
for a total of 8 visits using a total of 24 orbits (Fig. |2]i. A de- 
tailed summary of the observations is given in Table [1] For 
the assessment of the stellar populations in the general field 
of the LMC in the vicinity of the observed systems we used 
WFPC2 observations of the LMC field available in the HST 
Data Archive. We retrieved and analyzed images obtained 
in two campaigns, observed with parameters similar to ours. 
Details about these data are given also in Table [T] 

2.2. Photometry 

Photometry is performed using the package HSTp hot (Ver. 
1.1) specifically designed for WFPC2 imaging ("Dolphijil 
i2000.) . HSTphot is tailored to handle the undersampled nature 
of the point-spread function (PSF) in WFPC2 images and uses 
a self-consistent treatment of the charge transfer efficiency 
(CTE) and zero-point photometric calibrations. This pack- 
age has the ability to perform photometry simultaneously for 
all exposures in different filters. Another advantage of HST- 
phot is that it allows the use of PSFs which are computed di- 
rectly to reproduce the shape details of star images as obtained 
in the different areas of WFPC2 chips. For this reason, we 
adopt the PSF fitting option in the photometry routine, rather 
than use aperture photometry. Photometric calibrations and 
transformations were made according to Dolphin (2000), and 
CTE corrections were made according to Dolphin (200 2). We 
used the HSTphot subroutine mask to take advantage of the 
data quality files by removing bad columns and pixels, charge 
traps and saturated pixels. The same procedure is also able 
to properly solve the problem of the vignetted regions at the 
chip edges. The subroutine crmask was used for the removal 
of cosmic rays. 

The HSTphot photometry routine, hstphot, returns data 
quality parameters for each detected source, which can be 
used for the removal of spurious objects. After removing bad 
detections based on the quality parameters, we derived the fi- 
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Figure is omitted due to size limitations. 
Copy available upon request to the first author. 



Fig. 2. — Color composite images of the primary WFPC2 observations of our program GO-11547 constructed from the reduced images taken in the filters 
F555W (blue) and F814W (red). Top panel shows all Pointings 1, i.e., visits 1, 3, 5, and 7. Bottom panel shows the Pointings 2, i.e., visits 2, 4, 6, and 8. The 
corresponding cluster names are also given. 



TABLE 2 

Statistics on stars detected in filter-pairs with 
our photometry in each of pointings 1 . 



F300W F450W F555W F814W F656N 



LH 63 



F3()0W 


579 


565 


453 


515 


418 


F450W 




2 538 


2 389 


2 464 


500 


F555W 






5 027 


4 869 


394 


F814W 








5 753 


452 


F656N 










559 


LH60 


F3()0W 


348 


340 


294 


316 


257 


F450W 




2 509 


2 436 


2 467 


357 


F555W 






5 237 


5 050 


309 


F814W 








5 551 


332 


F656N 










409 


LH 72 


F3()0W 


445 


435 


362 


401 


282 


F450W 




2 022 


1 920 


1 963 


344 


F555W 






3 801 


3 698 


264 


F814W 








4 493 


306 


F656N 










421 


LH 88 


F3()0W 


206 


201 


173 


180 


141 


F450W 




1 777 


1 719 


1 737 


229 


F555W 






3 819 


3 715 


201 


F814W 








4 321 


207 


F656N 










285 



nal photometric catalogs of all stars found in at least two of 
the considered filters. The numbers of stellar detections with 
good photometric quality for each of the WFPC2 Pointings 1 
are given in Table|2] paired in all observed filters. From these 
numbers it can be seen that the richest stellar samples are de- 
rived in the filters F555W and F814W {V- and /-equivalent). 
This is due to the higher sensitivity of the camera in the cor- 
responding wavelengths in combination with the applied long 
exposure times in these filters. As far as Pointings 2 are con- 
cerned, observed only in these two filters, the numbers of stars 
detected in both filters with high photometric accuracy are 
4 397 in LH 63, 5 094 in LH 60, 3 886 in LH 72, and 3 958 in 
LH 88. 

We combined the photometric catalogs derived from the 
two pointings per cluster into one final catalog of stars 
detected in each targeted region. This process revealed 
12 511 stai-s in LH 60, 11 785 in LH 63, 9 756 in 
LH 72 and 10 597 in LH 88 detected in both F555W 
and F814W filter s. According to the methodology origi - 
nally presented by iGouliermis. Brandner. & HenningI (120061) . 
and subsequently applied to various star-forming regions in 
both of the Magellanic Clouds (see, e.g.. .Sabbi et al.. .20071 : 



Gouliermis et al. 2007; Schm alz et alJ l2008t ICignoni et aP 
2009i: iVallenari. Chiosi. & Sordoll2010h . deep HST imaging 
in V- and /-equivalent wavelengths provides an efficient de- 
tection of PMS stellar populations from their positions in 
the corresponding CMD. Therefore, we assess the PMS stars 
in the observed regions from the complete photometric cat- 
alogs of stars detected in both the F555W and F814W fil- 
ters in both pointings of each cluster The corresponding 
stellar charts of the regions are shown in Fig. |3] In addi- 
tion, since these filters provide the richest stellar samples, 
we also base the final multi-band photometric catalog of 
each cluster on the catalog of sources detected in both the 
F555W and F814W filters. A sample of this photometric 
catalog is shown in Table [3] Naturally, the photometry of 
stars identified in other pairs of filters are extremely useful 
and will be considered in subsequent analyses. For exam- 
ple the samples of stars detected in both F300W and F656N 
(U- and Ha-equivalent) certainly host PMS stars with accre- 
tion, which will be revealed from their U- and Ha-excess 
emission (see, e.g., iR omaniello. Robbe rto. & Pan agia 2004]; 
iDe Marchi. Panagia. & Romanielloi i2010l) ; an analysis that 
we will apply in a forthcoming study. Moreover, multi- 
wavelength coverage of bright main-sequence (MS) stars is 
essential for the assessment of the interstellar reddening, 
which we perform later in this study. 

2.3. Photometric accuracy and completeness 

Photometric uncertainty and incompleteness are two factors 
that strongly depend on the targeted region. As a consequence 
both are higher in the frames centered on the Hll nebulae, due 
to higher stellar crowding and confusion by diffuse nebular 
emission. Fig. |4]shows typical uncertainties of photometry as 
a function of the magnitude for all filters, as obtained from our 
photometry on Pointings 1. As seen in Fig. |4], photometric 
uncertainties behave similarly for each region in the F450W, 
F555W, and F814W filters. The F555W and F814W obser- 
vations reach much fainter magnitudes accurately than those 
with the F300W and F450W filters. On average, the uncer- 
tainties of our photometry have e < 0. 1 mag for TO555 ^ 25.4 
in all systems. 

The completeness of our photometry was evaluated on the 
basis of artificial star experiments performed with the native 
artificial star function built into HSTphot for the creation of 
simulated images by distributing artificial stars of known po- 
sitions and magnitudes. This utility allows the distribution of 
stars with similar colors and magnitudes as in the real CMD. 
It should be mentioned that due to the small image size of 
artificial stars generated within HSTphot the appearance of 
few bright saturated stars in each observed region leads to a 
somewhat underestimated completeness. The results of this 
process for the four broad-band filters used for Pointings 1 
of each cluster are shown in Fig. |5] From these plots it can 
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TABLE 3 

Sample of the multi-band photometric catalog of stars primarily detected in both F555W and F814W filters. 



Star 




R.A. 


Decl. 


X 


Y 


"1300 


£300 


"1450 


<:450 


"1555 


1:555 


m8i4 


£814 


"1656 


i:656 


# 




(deg 


J2000) 


(pixels) 






















1 


81 




—6/. 4254913 


149.25 


562.03 


20.489 


0.127 


19.280 


0.007 


18.407 


0.002 


17.205 


0.002 


17.361 


0.022 


2 


81 


959/015 


—6/. 4308624 


642.66 


148.52 


17.245 


0.008 


18.327 


0.003 


18.381 


0.002 


18.437 


0.003 


18.274 


0.030 


3 


81 


9557800 


—67.4288330 


466.22 


241.60 


99.999 


9.999 


20.181 


0.007 


19.108 


0.002 


17.397 


0.002 


17.737 


0.017 


4 


81 


9397507 


-67.4293671 


432.98 


729.03 


17.231 


0.009 


18.404 


0.004 


18.472 


0.002 


18.565 


0.003 


18.171 


0.033 


5 


81 


9583130 


-67.4300308 


571.80 


180.04 


17.182 


0.006 


18.321 


0.003 


18.402 


0.002 


18.477 


0.003 


18.214 


0.025 


6 


81 


9457169 


-67.4278641 


343.66 


531.86 


17.543 


0.008 


18.570 


0.003 


18.537 


0.002 


18.368 


0.003 


18.164 


0.026 


7 


81 


9394455 


-67.4292374 


421.46 


736.55 


17.805 


0.012 


18.637 


0.005 


18.623 


0.002 


18.643 


0.003 


18.251 


0.042 


8 


81 


9534912 


-67.4265442 


276.55 


282.40 


17.527 


0.009 


18.653 


0.003 


18.727 


0.002 


18.775 


0.003 


18.486 


0.028 


9 


81 


9536591 


-67.4270477 


316.44 


283.45 


18.109 


0.012 


18.827 


0.004 


18.852 


0.002 


18.818 


0.003 


18.609 


0.035 


10 


81 


9566727 


-67.4309464 


635.36 


240.53 


20.855 


0.079 


19.856 


0.007 


19.213 


0.003 


18.185 


0.002 


18.255 


0.024 


11 


81 


9431229 


-67.4300385 


501.62 


635.93 


18.301 


0.019 


19.129 


0.007 


19.098 


0.002 


19.022 


0.004 


18.708 


0.054 


12 


81 


9473572 


-67.4330444 


755.68 


545.32 


18.312 


0.019 


19.074 


0.005 


19.082 


0.003 


19.069 


0.004 


18.818 


0.056 


13 


81 


9619293 


-67.4312439 


682.93 


86.19 


17.757 


0.021 


18.759 


0.007 


18.822 


0.004 


18.882 


0.005 


18.577 


0.032 


14 


81 


9591370 


-67.4265823 


305.71 


113.45 


21.080 


0.072 


19.995 


0.013 


19.314 


0.004 


18.262 


0.003 


18.312 


0.024 


15 


81 


9490433 


-67.4320602 


687.11 


482.87 


18.534 


0.016 


19.196 


0.004 


19.176 


0.002 


19.176 


0.004 


18.942 


0.058 



Note. — In this example the 15 first records of the catalog for LH 63 is shown. '99.999 ± 9.999' values coiTespond to non-detection of the source in 
the specific filter. 
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Fig. 3. — Stellar charts of all stars detected with our WFPC2 photometry in both F555W and F814W filters in each of the obsereed Hll regions. Coordinates 
are given in seconds of arc from a reference point, which is selected to coincide with the center of WFP C2 Pointing 1 for every region. North is up and East is to 
the left of the charts. Bright candidate accreting PMS stars, identified from their Ha-excess (see ^ 14.11 . and OB-type stars found with our photometry in shorter 
wavelengths (§|4) are marked with red and blue symbols repsectively. 



Figure is omitted due to size limitations. 
Copy available upon request to the first author. 



Fig. 4. — Uncertainties of photometry as derived by HSTphot for all filters in the fields covered by Pointings 1 for each cluster. 

be seen that our photometry in Pointing 1 is more complete systematically for stars found in filters F555W and F814W. 
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Figure is omitted due to size limitations. 
Copy available upon request to the first author. 



Fig . 5 . — Completeness as it is evaluated by our artificial stai' experiments for the four broad-band filters in the regions covered by each Pointing 1 . 



In most of the clusters, the photometric completeness in the 
F300W filter drops at brighter magnitudes, msoo — 22, while 
in F450W it is somewhat fainter. On average, we reach the 
limit of 50% completeness at 771555 ~ 26. Completeness in 
both F555W and F814W in WFPC2 Pointings 2 is sUghtly 
improved due to lower crowding and confusion. 

3. INTERSTELLAR REDDENING 

In this section we quantify the effect of interstellar extinc- 
tion in the observed regions. We apply a statistical measure- 
ment of the visual extinction. Ay, in each observed region 
by constructing the distribution of the best-observed MS stars 
according to their extinction and obtaining its average value 
and 1(7 uncertainty. The calculation of Ay toward the first 
WFPC2 pointing of each system is directly achieved from our 
observations in all four broad-band filters through compar- 
ison of data points in color-color diagrams to a model MS 
in the WFPC2 photometric system (Girardi et al. 2002). For 
this calculation we considered only stars with S/N > 10 in 
all four bands, corresponding roughly to TO450 < 19.5 and 
— "1450 < 0. However, the second WFPC2 pointing of 
each system is observed only in the bands F555 W and F8 14W, 
and therefore we apply an additional independent Ay mea- 
surement based on the comparison of the observed CMD po- 
sitions of the upper MS stars to that expected according to the 
model MS. This measurement is applied for MS stars with 
S/N > 10 in both bands, corresponding to TO555 < 19.5 and 

771555 - "7,814 < 0.5. 

The extinction law depends on the nature and composi- 
tion of interstellar dust grains. Therefore the law will vary 
depending on environment such as Galactic star forming re- 
gions or for diff erent galaxies like the Large and Sma ll Mag- 
ellanic Clouds (iGochermann & Schmidt-Kaleill2002h . Sev- 
eral investigations of the reddening law in the LMC show 
that the average LMC extinction curve does not differ from 
that of the Milky Way at wavelengths covered by B-, 
V- and /-eq uivale n t filters (e.g ., Koo r nneef &_Cod^ 1981j 
Nandv et a l.] Il981t iFitzpatrickl fl9M. ISauvage & Vigrouxl 
~99k MisseltetalJll999l) . 3ut with differences appearing in 
the U band and for shorter wavelengths. Therefore, for the 
evaluation of extinction we assume the typical Galactic ex- 
tinction law (e.g., Cardelli et al. 1989; Fitzpatrick & Massa 
[T990I) . parameterized by a value of Ry = Ay /E{B -V) ^ 
3.1. Relative extinction was determined in the WFPC2 pass- 
bands from the analytical e xpression of stellar extinction of 
iFitzpatrick & Massal (119901) . 

Our results are shown in Fig. |6] where we plot the de- 
rived Ay distributions for each Pointing of every observed 
region. The average extinction values per region and pointing 
are given with their 1 a errors in Table |4] The fits applied to 
the histograms of Fig. |6] consider all stars within the selected 
S/N. It should be noted, however, that the CMD positions 
stars and most distant from the MS coincide with those of 
emission-line star-disk systems such as Be, B[e] stars (e.g.. 



TABLE 4 

Visual Extinction, toward the observed 

REGIONS. 



Region 


WFPC2 Pointing 1 


WFPC2 Pointing 2 


"UBVr 


"vr 


"V/" 


LH63 


0.49 ±0.19 


0.37 ± 0.09 


0.43 ±0.17 


LH60 


0.37 ±0.10 


0.31 ±0.12 


0.44 ± 0.21 


LH 72 


0.33 ±0.13 


0.29 ± 0.09 


0.41 ±0.15 


LH 88 


0.73 ± 0.64 


0.81 ±0.70 


0.76 ± 0.42 



Note. — The given ± values are the measured spreads, cr, 
in Ay, not the uncertainties in the measured values. 



IWisniewski et al.l2007|lKrausl2009l). orof Herbig Ae/Be stars 
(e.g.. iNishiyama et al.N2007l: IClayton et a l. 2010), which are 
expected in such star-forming regions. In addition, the ap- 
pearance of binaries among MS stars will naturally lead to a 
small overestimation of the evaluated reddening. As a conse- 
quence, the derived values and errors of Ay should be consid- 
ered as the upper limits for the actual extinction. It is worth 
noting that, as shown from the values of Table]?] while in gen- 
eral Pointing 2 is less crowded and less affected by the Hll 
emission than Pointing 1, Ay there is generally larger than in 
the latter This is due to the higher concentration of dust on 
the rim of the Hll region, as is typical in such star-forming 
regions, which produces higher extinction on the periphery of 
the nebula. 

In our subsequent analysis we will consider these upper 
values, rather than any lower average Ay derived from the 
distributions after trimming the higher values bins. We will 
specifically use the maximum Ay value, derived from the 
values of Table ]4] as (Ay)^^^^ + fJmax for each region. The 
reason for this approach is the fact that low-mass PMS stars 
in star-forming regions of the Magellanic Clouds occupy the 
fainter-redder part of the visual CMD, which c an be strongly 
affected by reddening (see, e.g.. iGouUermis et al. 2010). As 
a consequence faint low- and intermediate-mass MS stars 
may be confused for low-mass PMS stars due to reddening. 
Therefore, since it is important to apply a safe distinction 
between faint reddened MS and true PMS stars we quantify 
the maximum possible effect of visual extinction to the ob- 
served stellar samples in our qualitative distinction between 
PMS and reddened MS stars in the Hess diagrams and field- 
subtracted CMDs of the observed regions (see § ]5]l. On the 
other hand, the quantitative selection of PMS stars according 
to their cross-section CMD distributions is not affected by the 
considerations of maximum reddening values, because these 
distributions in the CMDs of the less affected regions LH 60, 
LH 63 and LH 72 peak at colors redder than those of the ex- 
pected most-reddened MS stars according to our reddening 
determinations (see §]6]). 

4. COLOR-MAGNITUDE DIAGRAMS 
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Fig. 6. — Visual extinction distributions as derived from upper MS stars witli S/N > 10 with the use of observations in all four broad-band filters in Pointings 
1 , and in the V- and /-equivalents for Pointings 2 for every region. We evaluate extinction in Pointings 1 also with the use of only V- and /-equivalent filters 
for reasons of completion and comparison to those for Pointings 2. The best-fitting distributions to the histograms are overlaid. The derived average (peak) Ay 
values and the corresponding la uncertainties are given in TableO [A color version of this figure will be available in the online journal.] 



The msoo — ?7T-450 versus 771450 (J/B-equivalent) and 771,555 — 
7ngi4 versus 777514 (^/-equivalent) CMDs of the stars de- 
tected in the areas of the clusters are plotted in Figs. |7]and 
[8] respectively. Since only Pointing 1 was observed in the U- 
and B- equivalent wavebands, the corresponding CMDs show 
the stellar populations included only in this pointing of each 
cluster. On the other hand, both Pointing 1 and 2 were ob- 
served in the V- and /- equivalent wavelengths, and therefore 
the corresponding CMDs show the stellar populations covered 
by both pointings. Reddening vectors corresponding to the 
highest Ay values derived in §[3] as they are given in TablelH 
are also plotted in the CMDs. In order to select stars with the 
overall best photometry, we have used the average error in the 
two bands, 62, considered for every CMD set, given as 



S2 = 



-300 



+ e: 



450 



S2 = 



-555 



■^814 



(1) 



respectively. The numbers of stars with 62 < 0.1 mag in every 
stellar catalog are given at the bottom-right of each CMD. 

In the CMDs of Figs. [7] and [8] it is shown that all observed 
regions are characterized by a variety of stellar species. The 
two sets of CMDs for every region, due to the selected fil- 
ters and exposure times, track different types of stars, which 
nevertheless belong to the same system. In particular, the 
UB CMDs of Fig. I2] indicate a sharp upper main-sequence 



(UMS) populated by the blue intermediate- and high-mass 
stars, which typically characterize embedded young LMC as- 
sociations. According the the ZAMS model (plotted over the 
CMD for LH 63) all observed clusters host well-populated 
UMS mass functions with the most massive stars being O- 
type dwarfs with masses between 20 and 30 Mq. These stars 
and their relation to the observed clusters will be the subject of 
another subsequent analysis of the massive stellar content of 
these clusters. The OB stars revealed in the observed clusters 
from the U B-equivalent CMDs are shown with blue symbols 
in the stellar charts of Fig.|3] 

In the present analysis we concentrate on the rich VI CMDs 
of the clusters (Fig. |8]l. These CMDs are populated by a 
larger variety of stellar sources, covering the MS (faint and 
bright), the turn-off (TO), the red clump (RC), the red giant 
branch (RGB) and the PMS. These features represent two sep- 
arate populations in the LMC, i.e., the evolved population of 
the general galaxy-field and the young and currently form- 
ing stellar populations of the clusters themselves. We de- 
scribe each of these features in the following text. It should 
be noted that the observations presented in Fig. |8] are mainly 
focused on the faint stars, and therefore long exposures were 
performed in these filters. As a consequence, the UMS and 
the tip of the RGB are saturated and thus not shown. All 
four CMDs of Fig. |8] share all aforementioned features, and 
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Fig . 7 . — The U _B-equivalent CMDs of the stars detected with our photometry in both F300W and F450W filters in Pointings 1 of the observed regions. Arrows 
represent the maximum Ay measurements given in Table|4] Typical photometric errors in both magnitudes and colors are shown on the left of each CMD. The 
numbers of stars with S2 < - 1 mag are given at the bottom-right of each diagram. The positions of these stars are indicated in the CMDs by thicker points. The 
ZAMS from the lGirardi et al.l J2002I) grid of models is overlaid with indicative stellar mass values for MS stars. [A color version of this figure will be available 
in the online journal.] 
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Fig. 8. — The y /-equivalent CMDs of the stars detected with our photometry in both F555W and F814W filters in both Pointing 1 and 2 of each observed 

region. The numbers of stars with 52 < 0.1 mag are given at the bottom-right of each diagram. The loci of the PMS stars cover the red sequence at m8i4 ^ 21, 
which is almost parallel to the MS. The intrinsic characteristics of these objects and possibly different ages among them are responsible for the spread of their 
CMD positions (see discussion in §[7] It is interesting to note that the CMD of LH 88 shows clear signs of extensive differential extinction in the system, apparent 
in the broad MS and elongated RC. 



specifically apart from the UMS there is a pronounced TO at 
mgi4 ~ 21 and a RGB with its RC located at nisu — 18 
and TO555 — msi4 — 1.0. These two features are known 
t o be typical of the LMC field, as we demonstrate later (in 
§ 14.21 see also Fig. [Tol i. Below the TO, moving to fainter 
magnitudes, the low main sequence (LMS) is the dominant 
feature of the CMDs. Previous HST imaging of other LMC 
star-forming regions has demonstrated that this population be- 
longs entirely to the general LMC field and not the stellar 



system s dGouUermis et alj l2007t IVallenari. Chiosi. & Sordol 
2010). On the other hand, the faint stellar membership of the 
young clusters at roughly the same brightness range with the 
LMS is located at redder colors in the CMDs, where a promi- 
nent broad sequence of PMS stars can be clearly seen in three 
of the observed regions. As we show later (§|5]) these stars be- 
long only to the star-forming regions, and they are not related 
to the general LMC field. 
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Fig, 9. — Left: Diagnostic for stars with Ho excess emission in the _R — Ho, V — /-equivalent color-color diagram. In this diagram the positions of normal 
stars without Hq emission, i.e., field or weak-Une T Tauri stars, are indicated by the blue line, and reddened field stars by the green. The group of stars with 
R — Hq > 1.0 mag are mostly stars, which exhibit Ha excess (plotted with red symbols). Right. The VZ-equivalent CMD of the whole sample of observed stars 
in all four regions with Ho excess emission candidate stars overlaid with red thick symbols. In this plot the CMD positions of bright PMS stars with Ho excess 
can be readily identified, but the faint PMS populations are not covered due to the low sensitivity of the Ho filter The nature of the detected bright blue sources, 
as well as of cool objects with Ho excess will be investigated in a subsequent study. [A color version of this figure will be available in the online journal.] 



4.1. Stars with Ha excess 

PMS stars exhibit excess in their Ha emission due to ac- 
cretion. As a consequence the appearance of Ha emission 
from stars in the clusters is a direct evidence of the exis- 
tence of PMS stars in them. Nevertheless, if the age 
of the observed clusters is of the order of 3 M yr, Ha 
emiss ion stars should be very rare (see, e.g.. Su ng et al.l 
12004ft . In addition, the accretion rate of most PMS stars 
in the clusters may be very low, and therefore only few ac- 
tive stars would show Ha emission, because only stars with 
equivalent widths W^(Ha) > 10 A can be detected photo- 
metrically with confidence jSung et al.ll2008l) . Sources with 
Ha excess emission can be photometrically identified by the 
comparison of th eir Ha and i?-equivalent magnitudes (e.g., 
ISung etal.lfT99^ . As our dataset lacks observations in the 
i?-equivalent WFPC2 filter F675W, we determined a "syn- 
thetic" mgys magnitude for every star derived by interpolation 
be tween data in the F555W and F8 14W bands, as described 
bv lDe Marchi. Panagia. & Romanie llo (2010). In total, 2 263 
stars are detected in all three F555W, F814W and F656N 
bands with the best photometric quality in all four observed 
regions. We identify the Ha excess stars in our sample by 
plotting the mgrs — mgse color index against the 771555 ^ ^7t-814 
(Fig. |9] left). This diagnostic is equivalent to the use of the 
mean magnitud e of V and / as a pseudoc ontinuum magnitude 
at Ha (see, e. g.. ISung et a"ni2000l: ISunF & Bessell 2004). 

Considering that the reddening vector in Fig. |9] (left) is 
nearly parallel to the 771555 ~ "1814 axis, the clump of stars 
around 771555 ~ ^814 ^1.0 and 771575 ~ W656 ~ 0.3 rep- 
resents highly reddened background stars. The green line in 
the fi gure indicates the meanline of these non-em ission stars 
(e.g.. lPanagia"etani2000l: iSung & Bessell l2004h . In Fig. |9] 



(left) a large scatter in the 777,675 — "1,656 index of normal 
stars, i.e., field stars or weak-line T Tauri stars, can also be 
seen. Taking into account the positions of these stars, indi- 
cated in the figure by the blue line, it is safe to select stars 
with 777675 ~ "7,656 ^ 1.0 as candidate Ha emission stars. 
As a consequence, for a first-order identification of the best 
candidates we apply tentatively the selection criterion requir- 
ing all Ha emission stars to have 777575 — fnese > 1-0 mag. 
Objects bluer th an 777555 — 7?78i4 — 0-2 are p robably Ae/Be 
stars ( De Marchi. Panagia. & Romanielloll2010i) . but consid- 
ering the youthfulness of the clusters we should also expect 
a few massive PMS stars, i.e., Herbig Ae/Be stars. This 
selection revealed 78 sources, which appear in red in the 
"1-675 ~ 7r7656, 777555 — "7814 diagram. These sources are 
marked with red symbols in the T^/-equivalent CMD of all 
detected sources shown in Fig. |9] (right), and in the stellar 
maps shown in Fig. [3] 

While our F656N observations provide direct evidence that 
the observed clusters host accreting PMS stars of intermedi- 
ate and high mass, they are designed for the investigation of 
the bright stellar content of the clusters, and therefore they 
are not deep enough to cover the more populous low-mass 
PMS stellar content of the clusters. In fact, the aforemen- 
tioned process identified PMS accreting candidates down to 
only 7778 1 4 ~ 22, much brighter than the vast majority of the 
low-mass PMS stars in the clusters. For the statistical iden- 
tification of the latter we make use of the rich ^/-equivalent 
CMDs, as described in the following sections. Nevertheless, 
a further thorough investigation of bright PMS stars with ac- 
cretion, based on their Ha equivalent widths, as well as their 
excess emission in the J7-equivalent filter, and the analysis of 
the relation between Ha strength and U -excess with accretion 
rate will be presented in a forthcoming paper. 
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4.2. Contribution of the general LMC Field 

Naturally, the CMDs of Fig. |8] are 'contaminated' by the 
stellar populations of the general field of the LMC. In order to 
assess the field stellar population in the area of the observed 
regions, and evaluate this contamination, we retrieved from 
the HST Data Archive WFPC2 observations of empty regions 
of the LMC. Special care was taken to select observations, 
which are the most representative of the LMC field in the gen- 
eral vicinity of our LMC 4 target regions. In addition the fields 
should preferably have negligible differential extinction, and 
observational setup comparable to our own. Two WFPC2 
pointings covered in the F555W and F814W bands matched 
best these criteria. The fields are observed within HST pro- 
grams GO-10583 (PI: C. Stubbs) and GO-10903 (PI: A. Rest), 
both focused on the detection of LMC microlensing events in 
Cycles 14 and 15 respectively. Details on these observations 
are given in Table [T] Both our clusters and the general field 
were observed with almost identical setup of the instrument 
and exposure times, allowing a direct comparison between the 
VI CMDs of the clusters and the local LMC field. 

The CMDs of these LMC control-fields, shown in Fig.fTOl 
outline the complete star formation history (SFH) of the 
galaxy, which has been well constrained by earlier studies 
with WFPC2 imaging. The results of these studies provide 
an accurate account of the evolutionary status of the stars 
we see in the CMDs of Fig. [TO] According to these results, 
in general, the LMC disk is characteri zed by a continuous 
SFH for the last ^ 10 to 15 Gyr (e.g., iSmecker-Hane et al.l 
120021) . with events of enhanced star formation that 
took place betwee n 1 and 4 Gyr ago (Gallagher et al.' 
1996; Elson. Gilmo re. & Santia go 1997; Gehaetal. 1998; 
Castro et al.l 120011: iJaviel. Santiago. & KerbeH 120051) . These 
events are well recorded in the control-field CMDs of this 
study, as shown by the isochrones overlaid on the overall field 
CMD of Fig, nn (Jeft). These isochrones correspond to ages 
of roughly 0.8, 1.25, 2, 3, and 8 Gyr, proving this area to 
be somewhat younger than those described in the aforemen- 
tioned studies, which are located outside LMC 4. A negligi- 
ble visual extinction of Ay — 0.075 mag, and the metallicity 
(Z = 0.008) and distance of the LMC (to - M = 18.5 mag) 
are applied to these models. 

5. IDENTIFICATION OF PMS STARS IN THE OBSERVED REGIONS 

Low-mass PMS stars in LMC star-forming regions are 
known to cover the red-faint part of their ^/-equivalent 
CMDs (see, e.g., Vallenari, Chiosi, & Sordo 2010), and the 
regions observed here do not seem to be an exception. In- 
deed, a visual comparison between the CMDs of Fig. |8] and 
those of Fig. [TO] verifies the existence of a broad sequence 
of PMS stars almost parallel to the LMS in the CMDs of 
the star-forming regions, which is completely missing from 
the control-field CMDs. In addition, the reddening range 
of early-type stars in the star-forming regions (Table 2}, in 
comparison to that of the control-field {Ay — 0.075) shows 
that reddening alone cannot explain the existence of very red 
population almost parallel to the LMS in the CMDs of the 
observed regions. In this part of the analysis we focus on 
methods for the identification of these PMS in the observed 
regions. In the cases of LH 60, LH 63, and LH 72 the appear- 
ance of PMS stars in the CMD is rather clear, but in the LH 88 
region the higher reddening produces confusion in the unam- 
biguous identification of PMS stars. The radical difference 
in the comprised populations between clusters and field is 
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Fig. 10. — The y/-equivalent CMDs of the stars detected with our pho- 
tometry in both F555W and F8I4W filters in the two selected empty areas for 
the assessment of the general field of LMC in the region of LMC 4. These 
WFPC2 images were retrieved from the HST Data Archive (see §|2]and Ta- 
ble[T). The numbers of stars with 82 < 0.1 mag are given at the bottom-right 
of each diagram. 

demonstrated more clearly in Fig. [TT] where next to the over- 
all control-field CMD with overlaid isochrones for evolved 
populations (left panel), the CMD of LH 60 is plotted with 
overlaid isochrones for PMS populations (right panel). These 
PMS models are constructed with t he Frascati Raphson New- 
ton Evolutionary Cod e (FRANEC. I Chieffi & StranTerollT989t 
iDeglTnnocenti et al.l2 008) for the metalhcity of the LMC and 
the WFPC2 photometric system. They correspond to ages of 
0.5, 1, 2, 3, 4, 5, 7.5, 10, 15, 20 and 30 Myr This broad age- 
coverage demonstrates a well-documented problem in such 
observations, i.e., the wide spread of faint PMS stars, which 
does not necessarily correspond to any real age-spread among 
these stars, but it may well be the result of biases introduced 
by observational constrains, i.e., photometric accuracy and 
confusion, and the physical characteristics of these stars, such 
as variability, binarity and circumstell ar extinction (see, e.g., 
iDa Rio. Gouliermis. & Gennaroll2010i for a detailed discus- 
sion). 

The visual comparison of the CMDs of Fig. [TT] outline 
three important aspects concerning the observed stellar popu- 
lations: ( 1 ) The completely different evolutionary status of the 
field stars from that of the stellar content of the young clusters. 

(2) The great difficulty in distinguishing the evolved from the 
PMS populations at magnitudes brighter than msu — 21.5. 

(3) At fainter magnitudes the heavy contamination of the clus- 
ter CMD by the field populations, particularly in the LMS. 
Taking into account these aspects, it is important to establish 
an effective methodology for the determination of the true 
PMS stellar population of each observed young cluster/star- 
forming region. In this section we present such a scheme 
with the application of a qualitative approach based on the 
use of differential Hess Diagrams, and of a statistical field- 
subtraction technique based on the Monte Carlo method. Al- 
though these techniques provide a first distinction of the red- 
der PMS from the bluer MS low-mass population of the clus- 
ters, they are indirect approaches in the determination of the 
true PMS populations. In particular, with each of these meth- 
ods one can assess the field MS stellar contribution in the ob- 



Pre-Main-Sequence Stars across Shapley Constellation 111 in the LMC 



11 




-1 1 I J -t 1 2 3 
F555W-FB14W F555W-FBi4W 



Fig. 1 1 . — Left: VZ-equivalent CMD of the stars in the general LMC field 
at the vicinity of LMC 4 with isochrones for evolved stars overlaid. The mod- 
els are taken from the Padova grid (Girardi et al. 2002) for the LMC metal- 
hcity (Z = 0.008) and correspond to ages between roughly 1 and 10 Gyr 
in accordance with the established SFH of the galaxy. Right VZ-equivalent 
CMD of the stars in the observed region of LH 60 with isochrones for PMS 
stars overlaid. The models are constructed with the FRANEC code in the 
WFPC2 photometric system and the metallicity of the LMC (see § [5). [A 
color version of this figure will be available in the onhne journal.] 

served CMD in terms of stellar densities in the CMD, rather 
than on a star-by-star basis. Therefore, with these methods 
we cannot define the actual stars that should be considered as 
the most probable PMS stars, except for the reddest ones. As 
a consequence, we utilize these methods in order to demon- 
strate in a qualitative manner the existence of low-mass PMS 
stars in the observed regions. A quantitative statistical de- 
termination of the actual membership of each PMS candi- 
date will take place in the next section with the construction 
of the stellar distributions of the observed stars along cross- 
sections of the faint part of each CMD. With this method a 
PMS membership probability is assigned to each red star, ac- 
cording to its CMD position with respect to the cross-sections 
stellar distributions across the CMD. In our analysis we con- 
sider only stars with the best photometric measurements, i.e., 
with 62 < 0.1. 

5.1. Differential Hess Diagrams 

Fig.[T2]shows the ^/-equivalent Hess diagrams of the stars 
in the observed regions, as well as that of the overall gen- 
eral LMC field, i.e., the combined stellar content of the two 
control-fields. These diagrams are produced by binning the 
corresponding CMDs in 0.15 mag intervals in 7/1555 — ?7i8i4 
and 0.17 mag in TO814- The size of the binning is selected 
to be small enough to reveal the fine-structure in the stellar 
density, but also large enough to construct a realistic repre- 
sentation of the observed CMDs. Very small bin sizes would 
produce 'noisy' Hess diagrams, i.e., high numbers of small 
density peaks distributed throughout the diagrams, while large 
bin sizes would produce unrealistically thick Hess diagrams, 
which would smooth out any apparent effect of reddening or 
any true positional spread of stars in the CMD. The result- 
ing two-dimensional distributions were smoothed by inter- 
polation with a boxcar average. From the Hess diagrams of 
Fig. [12] it can be seen that the highest concentration of stars 
in all observed regions appears in the LMS. This part of the 



CMD i s we ll-known to belong entirely to the general LMC 
field (§11211. In the cases of LH 60, LH 63 and LH 72, at red- 
der colors from the LMS and almost parallel to it, there is an 
apparent lower concentration of stars, which does not appear 
in the field. This feature corresponds to the faint PMS stars of 
the star-forming regions. 

The differences between the observed systems and the field 
are more apparent in the residual or differential Hess dia- 
grams. These are constructed with the subtraction of the field 
Hess diagram from those of the observed regions. Before 
this subtraction, we simulated the CMD of the observed re- 
gions by shifting fraction of stars in the CMD of the control 
field according to the reddening law, so that the morphol- 
ogy of the RGB in the simulated CMD is similar to that of 
the observed ones. Then we adjusted the number of RGB 
stars of the simulated field CMD to that of the target region. 
In addition, a broadening of the MS band due to photomet- 
ric error was also simulated using a Monte Carlo technique. 
The resultant differential Hess diagrams are shown in Fig. [13] 
These diagrams for LH 60, LH 63 and LH 72 highlight the 
characteristic features of the V^/-equivalent CMDs of typi- 
cal LMC star-forming regions, i.e., a prominent UMS down 
to the intermediate-mass regime and a broad sequence of red 
faint PMS stars reaching the sub-solar mass-regime. In the 
differential Hess diagrams of Fig. [T3] there are some resid- 
ual stellar concentrations remaining in particular at the RC 
and some locations of the RGB. They are most likely due to 
the small number of stars in the corresponding grids, which 
do not allow a proper subtraction from the CMD of the ob- 
served regions. In addition, in the differential Hess diagrams 
of Fig. [13] there are no intermediate-mass PMS found in the 
Kelvin-Helmholtz contraction phase, located between the up- 
per PMS and the MS, although there should be some such 
stars at least in LH 60, LH 63 and LH 72. A possible reason 
for this lack is the over subtraction of field stars on the TO. 
These points suggest that the differential Hess diagrams are 
subject to the normalization of the stellar numbers of the field, 
which nevertheless does not represent the real background 
field in the observed regions, as well as to the simulation pro- 
cess for the reproduction of this real field. Therefore, we rely 
on the Hess diagrams only for demonstrating the existence of 
a true PMS population of each region, and not for quantifying 
its members. 

Considering the stellar distributions in the red-faint part of 
the CMDs of LH 60, LH 63 and LH 72, while the appearance 
of a sequence of stars in this part of the CMDs is quite obvi- 
ous (see Figs [8j [12] and [T3] l, one may argue that this part of 
the CMD is populated as a result of the effect of differential 
reddening on MS stars, rather than by actual PMS stars. How- 
ever, if we consider the visual extinction measurements per- 
formed in §[3]and the maximum of the derived values given in 
Table [4] as well as the corresponding uncertainties, we can- 
not fully assign the observed red-faint stellar sequences to 
reddening of the MS alone. In particular, while the bright- 
est part of the red-faint stellar sequence, at nisu — 21 and 
"1555 — mgi4 ~ 1, may be partly explained by extinction, 
the colors covered by the faint part of the sequence are com- 
pletely independent from a hypothetical differential reddening 
of the MS, as well as from the measured photometric errors at 
these magnitudes. This can be directly verified on the CMDs 
shown in Fig.[8]by considering the directions of the reddening 
vectors. 

The analysis thus far shows that the PMS feature in the 
CMDs and Hess diagrams of the observed clusters appears 
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Fig. 12. — Hess diagrams of the observed regions, and the general LMC field, constracted from their ^/-equivalent CMDs. The diagrams are color-coded 
according to their stellar number densities indicated in the color bars. [A color version of this figure will be available in the online journal.] 



to be quite distinct from a well-defined MS in the cases of 
LH 60, LH 63 and LH 72. On the other hand, the case of 
LH 88 seems to be quite different, because 1) we do not ob- 
serve the same sequence of stars in the red faint part of the 
observed CMD and the corresponding Hess diagram of this 
region, and 2) the MS of LH 88 as seen in these diagrams 
is quite broader than those of the other regions and the field. 
This behavior can be explained by the higher extinction and 
its differential nature in the vicinity of LH 88. Indeed, there 
is a dense, thin residual sequence of stars to the red part of 
the LMS in the differential Hess diagram of LH 88 (Fig.fTsTl. 
which does match the loose broad sequence of PMS stars seen 
in the other regions. This sequence could be misinterpreted as 
an 'old' PMS population, because the extinction measured for 
LH 88 (Table|4l) can easily span its width. This stellar feature, 
thus, may include both reddened LMS stars and PMS stars. 
While we can not exclude the presence of PMS stars in LH 88, 



their loci should be well contaminated by reddened MS stars, 
and therefore their identification is not as straightforward as 
in the other regions. 

5.2. Statistical Identification of PMS stars: Subtraction of the 
field stars contamination 

The stellar populations observed in the regions of the clus- 
ters are subject to strong contamination by the general LMC 
field. This contamination can be assessed statistically with 
the use of CMDs of the LMC field (Figs.[Tni[n]Jeft) and sub- 
tracted from the cluster CMDs with the application of a Monte 
Carlo method. For this method we consider a set of elliptical 
subregions around each star in the CMD of each cluster and 
the same set of subregions in the CMD of the LMC field. We 
then statistically subtract from the CMD of each system the 
corresponding number of randomly selected stars in the field 
CMD for each of these subregions. We perform 100 itera- 
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Fig. 13.— Differential Hess diagrams of the observed regions derived by the subtraction of the Hess diagram of the general LMC field from those of the 
observed regions. [A color version of this figure will be available in the online journal.] 
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Fig. 14. — Field-subtracted y/-equivalent CMDs for the observed star-forming regions. The remaining stellar members are color -coded according to their 
membership probabilities derived from the repeated application of our Monte Carlo method for the field decontamination (see §|5]2). [A color version of this 
figure will be available in the online journal.] 



tions of the Monte Carlo subtraction for each cluster, and we 
assign, thus, a percentage of membership probability, p, to 
every remaining star-member The resultant field-subtracted 
CMDs of the clusters are shown in Fig. [14] overlaid on the 
originally observed CMDs (in grey). The stars remaining af- 
ter the subtraction process are color-coded according to their 
assigned membership probability with p > 55%. In this fig- 
ure it is shown that, as was the case in the Hess diagrams of 
the previous section, stars located in both the faint red and 
bright blue part of the CMDs are flagged as true members of 
the clusters. As it is observed also in the Hess diagrams of 
Fig. [13] after the field-subtraction there are few red evolved 
stars, which are falsely flagged as cluster members. These are 
residuals of the subtraction technique, which was not possi- 
ble to be removed due to natural inconsistencies between the 
CMDs of the regions and that of the control-field. 

The large majority of residual cluster-member stars with 
membership probability p > 95%, occupy the PMS part of the 



CMDs in Fig. [14] Naturally, this behavior characterizes the 
cases of LH 60, LH 63 and LH 72, where the PMS part of the 
CMD is fully populated by cluster-member stars down to the 
faintest detected magnitudes. On the other hand, in the case 
of LH 88, residual red stars with high membership probabili- 
ties correspond to somewhat brighter magnitudes and form a 
sequence connected to the residual field stars of the sub-giant 
branch and RC. Considering the high extinction of this region, 
the appearance of these red sources in a thin sequence almost 
attached to the red part of the LMS can be best explained by 
the effect of differential reddening on MS stars, rather than by 
these stars being PMS stars. The case forLH 63, which is also 
affected by - nonetheless lower - extinction, is not similar to 
that of LH 88, because the stars with high-membership prob- 
ability populate completely the PMS part of the CMD as in 
the cases of LH 60 and LH 72. It is interesting to note that in 
the cases of these two clusters, which show the lowest extinc- 
tion, there is a clear gap between the LMS and the PMS stars. 
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Fig. 15. — Example of the stellai' counting process for the construction 
of the distribution of stars along cross-sections through the observed VI- 
equivalent CMD of LH 60. The six selected cross-sections perpendicular to 
the MS and the PMS locus are indicated by the red lines, covering the whole 
faint part of the observed magnitude range. The strips selected for the binning 
of the stars along the cross-sections, and the construction of the correspond- 
ing stellar distributions are indicated by the blue lines, almost parallel to the 
MS. [A color version of this figure will be available in the online journal.] 

making the distinction of the latter more straightforward. 

Nevertheless, it should be noted that the results of the 
Monte Carlo subtraction technique are subject to its param- 
eters. In particular the resulted 'clean' CMDs and the remain- 
ing stellar numbers depend on the size of the individual el- 
liptical regions around each star considered for the random 
elimination of the field stars. Also this technique makes use 
of a general control-field rather than the actual field at the area 
of each observed region. Therefore, we treat the results of the 
Monte Carlo technique qualitatively and not quantitatively, as 
we did for the Hess diagrams. A quantitative analysis for 
the determination of the most probable PMS member-stars of 
each region is performed in the next section. Taking into ac- 
count the results of this section, and in order to eliminate the 
appearance of reddened field stars in our PMS samples, in our 
subsequent analysis we concentrate on the fainter part of the 
CMDs with 771814 > 21. 

6. MEMBERSHIP DETERMINATION OF PMS STARS IN THE 
OBSERVED REGIONS 

In this section we perform a statistical determination of the 
membership of the PMS stars in each observed star-forming 
region by constructing the stellar distributions along cross- 
sections of the observed CMDs. Since, as discussed above, 
at magnitudes brighter than msi4 ^ 21 there is a confusion 
between PMS and evolved stars, we limit our treatment to the 
fainter part of each CMD. We consider again only the stars 
found with photometric uncertainties S2 < 0.1. 



6.1. Stellar distributions along cross-sections of the CMD 

A more quantitative analysis of the observed PMS pop- 
ulations is achieved by counting the number of stars as a 
function of their (7^,555 — nisii) colors along a series of 
CMD cross-sections perpendicular to the PMS locus. This 
method is previously established for the identification of PMS 
stars in VI CM Ds of G alactic star-forming regions (see, e.g., 
ISherry. Walter. & WoUd 1200 4). An example of the counting 
process for the construction of the stellar distributions along 
six cross-sections, selected to cover the magnitude range of 
interest, i.e., 21.5 ^ 7^814 ^24.5, is shown in Fig. [TS] 
for LH 60. The cross-sections are indicated by the red lines 
with their numbers marked. Stars are binned in strips per- 
pendicular to the cross-sections and thus almost parallel to 
the MS. These strips, which are also almost parallel to typi- 
cal P MS isochro nes derived fr om PMS evolutionary model s 
(e.g- lPallafeSTa hler 1999; Sie ss. Dufour. & Forestinill2000l) . 
are indicated in Fig.[T5]by the blue lines. 

The constructed stellar distributions are shown in Fig. [16] 
Stellar numbers are corrected for incompleteness according 
to our completeness measurements described in § 12.31 Since 
the distribution of stars and nebula is not uniform across each 
observed field, confusion may also vary across each region, 
being higher at the most crowded regions. As a consequence, 
photometric completeness not only is a function of the ob- 
served magnitudes but also depends on the position of each 
star across the field. Therefore, we corrected the numbers of 
stars per bin in every strip on a single-star basis according 
to the magnitude and position of each star in the observed 
field. In Fig. [16] it can be seen that the distribution of stars 
for each cross-section is peaked on the field MS, which com- 
prises the large majority of the observed stars in each target. 
Our aim is to identify the stellar distributions that correspond 
to the redder PMS population of the star-forming regions. In- 
deed, as it can be seen for the regions LH 60, LH 63 and 
LH 72, their stellar distributions are not symmetrically peaked 
on the MS and are clearly extended to the red, especially for 
cross-sections with ?Ti8i4 ^ 24.0. It should be noted that, the 
red-ward asymmetry of their stellar distributions cannot be 
attributed to reddening alone based on our visual extinction 
estimates. 

For comparison to these distributions we show in Fig. [17] 
the corresponding stellar distributions derived from the CMD 
of the control fields (Figs. [10] fTTT i. From the distributions of 
Fig. [it] it is shown that field stars, in contrast to those of the 
star-forming regions, are distributed symmetrically around the 
MS. Again, the case of LH 88 fits more this behavior, rather 
than that of the other three regions of our sample. The stellar 
distributions in LH 88, shown in Fig. [16] appear symmetri- 
cally peaked on the MS, as in the case of the LMC field, but 
much broader than the general field, most probably due to 
the high extinction that characterizes this region. Under these 
circumstances, even if there are PMS stars in LH 88, we can- 
not accurately assess their numbers. In all of the discussed 
distributions, the wider spread of stars for cross-sections cor- 
responding to fainter magnitudes can be naturally attributed 
to the larger photometric uncertainties. 

Considering the above, the distribution of stars along cross- 
sections through the CMDs of the star-forming regions can be 
well represented by the sum of two distributions; one for the 
field MS stars and one of the native PMS stars of the regions. 
For demonstration we apply such a fit to the cumulative stellar 
distribution for each region and the field. These distributions 
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Fig. 16. — Stellar distributions along selected cross-sections of the CMDs of the observed regions. Since we are interested in the faint PMS populations of 
the clusters we constrain our analysis to six cross-sections through the CMD, covering the magnitude range 21.5 S; mgn $ 24.5. Stars are counted in bins 
0.035 mag wide and the derived distributions are plotted with respect to the coiTesponding average color (mrjss ~ "1814) of Each bin. Different colors are used 
for the distributions of different cross-sections. Stellar numbers are corrected for incompleteness. These plots reveal the PMS populations of the observed clusters 
as a secondary distribution red-ward from that peaked on the field MS stars. [A color version of this figure will be available in the onhne journal.] 
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In this manner, the peak of the distribution will keep following 
the original peaks on the MS of the individual cross-sections. 
On the other hand, since the counting bins (strips) are not ver- 
tical to the color axis, each bin in the cumulative distribution 
does not correspond to a single color value. Therefore, we as- 
sign to each bin the average color index, which roughly cor- 
responds to that of the average magnitude of msi4 — 23. 

We fit the number of stars per bin of the cumulative distri- 
butions to a double Gaussian function of the form 
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Fig. 17. — Same as Fig. ll6l for the total population of the general LMC 
field. [A color version of this figure will be available in the online journal.] 

are constructed by counting the stars in the selected strips par- 
allel to the MS for the whole considered magnitude range, i.e., 
in all six selected cross-sections, as they are shown in Fig.fTSl 



The first term describes the distribution of field MS stars and 
the second that of the PMS stars. We use a least-squares mul- 
tiple Gaussian fit performed by the interactive IDL routine 
XGAUSSFIT (by Don Lindler) to solve for six parameters, Ni, 
fii, and (Ti, where i = 1 is for the field and i ~ 2 for the 
PMS stars. The constructed cumulative stellar distributions 
across the CMDs of the observed regions and the LMC field 
are shown in Fig. [18] The best-fitted Gaussian distributions 
derived from our fitting process are also drawn; the total fit 
is indicated with a red line, while each of the two Gaussian 
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Fig. 18. — Cumulative stellar distributions along all considered cross-sections of the CMDs for the magnitude range 21.5 ^ mgi4 ^ 24.5. Stars are counted 
in bins 0.035 mag wide and their numbers are corrected for incompleteness. The distributions are plotted with respect to the average color (msss — msi4) 
of each strip (bin). The best-fit double Gaussian function is overlaid with a red line. Each component of this function is also drawn with orange lines. EiTors 
represent Poisson statistics. Residuals are plotted below each fit, demonstrating the goodness-of-fit. [A color version of this figure will be available in the online 
journal.] 



components is drawn with an orange line. The correspond- 
ing residuals of the fits are accordingly drawn in order to show 
the goodness of each fit. 

From the plots of this figure for LH 60, LH 63 and LH 72 
one can see that the red component of the stellar distributions 
along the cross-sections of the CMDs is indeed important and 
successfully represented by a Gaussian distribution. In the 
case of LH 88 the contribution of the red component to the 
stellar distribution turns out - as expected - not to be signif- 
icant, as can be seen by the fits. The cumulative distribution 
of stars in the field CMD is shown also in Fig.[T8]for compar- 
ison. Naturally, only a single-Gaussian fit could be applied to 



this distribution. The coefficients of the fits shown in Fig. [18] 
are given in Table |5] The values of this table show the re- 
markable coincidence of the primary peak in the distributions 
of the star-forming regions with the Gaussian function rep- 
resenting the general field. This coincidence, which is more 
prominent in LH 60, LH 63 and LH 72, provides additional 
evidence of the field origin of the MS population observed in 
the CMDs of these star-forming regions. 

This fitting process is performed for the stellar distributions 
along each of the selected cross-sections for each observed 
CMD. The distribution of the field stars per se is not interest- 
ing, but it is crucial because the number of PMS stars depends 
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TABLE 5 

Coefficients of the best-fitting double Gaussian 
function to the cumulative distributions of faint 
stars across the cmds. 





Gaussian component 1 


Gaussian component 2 


Region 
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^i2 FWHM2 


LH 63 


547 


0.82 


0.27 


162 


1.08 0.68 


LH 60 


888 


0.80 


0.24 


118 


1.03 0.78 


LH 72 


663 


0.81 


0.23 


96 


1.08 0.94 


LH88 


483 


0.89 


0.36 


93 


1.03 0.71 


Field 


1071 


0.82 


0.27 







upon the extrapolation of the field star distribution through 
the PMS locus. The total number of PMS stars along a cross- 
section is determined by the normalization of the Gaussian 
fit to the PMS distribution, N2,j (where j = 1, ...,6 is the 
cross-section number), and by (T2.j, which defines the width 
of the Gaussian. The full width at half maximum (FWHM) 
of the PMS population is 2(72 j \/2ln{2). The fi2.j parameter 
is the color 771555 ^ "1814 of the peak of the PMS distribution 
along a cross-section. Since each cross-section is a predefined 
line across the CMD, /i2.j specifies the magnitude that corre- 
sponds to the peak of the PMS distribution. The coefficients 
of the best-fit second Gaussian component, representative of 
the PMS population, for each cross-section are given in Ta- 
ble |6] 

We fit the field and PMS distributions jointly, and we de- 
rive the total number of PMS stars from the second Gaussian 
component of each fit by assigning a membership probabil- 
ity to each star. This probability is determined by the ratio 
of the second Gaussian component of the fit to the total fit 
along each cross-section through the CMD, i.e., the sum of 
both Gaussian components of the fit. Membership probabil- 
ities calculated in this manner vary from 0% to the blue of 
the PMS locus to >90% near the peak of the PMS distribu- 
tion and for redder colors in all cross-sections. We thus derive 
1 815 stars in LH 63, 1 230 in LH 60, 1 223 in LH 72 and 
only 314 in LH 88 with probabilities of cluster membership 
p > 95%. These numbers are completeness corrected. 

This particular statistical technique for the determination of 
cluster membership probabilities of PMS stars is quite self- 
consistent, because it makes use of the LMC field as it is ob- 
served within each region and not of a generic remote field 
as i n the case of the Monte Carlo field subtraction presented 
in § 15.21 Moreover, the use of the cross-sections distributions 
for the determination of the most probable PMS stars is inde- 
pendent of the reddening measurements in each region. It is 
interesting to note that this technique returned two times more 
PMS member candidates (with p > 95%) than the Monte 
Carlo field subtraction, except in the case of LH 88, where the 
derived numbers are comparable. The reason for this differ- 
ence lies on the elimination by the Monte Carlo technique of 
the red wing of the LMS from the CMDs of the star-forming 
regions as field candidates. Since this technique is based on 
the iterative random subtraction of candidate field contami- 
nants within individual CMD-regions around each observed 
star, the subtraction of different stars in each repetition results 
in low membership probabilities for all stars in CMD-regions 
where field stars are expected to be located. On the other 
hand, the cross-sections technique, being dependent solely on 
the actual CMD positions of the observed stars and their dis- 
tributions, returns a purely probabilistic membership determi- 
nation for each star, independent of the appearance of field 
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Fig. 19. — Cumulative stellar distributions across the CMDs of the ob- 
served regions for the stars found with membership probability p > 95%. 
The PMS populations of LH 60, LH 63 and LH 72 demonstrate an extraor- 
dinary similarity, which possibly indicates their common characteristics. On 
the other hand the distribution for LH 88 demonstrates the difficulty in disen- 
tangling the true stellar population of the region, since stars with high mem- 
bership probability seem to represent reddened MS stars and not PMS stars, 
which are possibly represented by the second (faint) peak. [A color version 
of this figure will be available in the online journal.] 

Stars or not in its surrounding CMD-region. 

7. DISCUSSION 

In Fig. [19] we show the cumulative distributions along the 
CMD only for the PMS stars of each region, i.e., for those 
stars found with a membership probability p > 95%. These 
distributions are constructed as those for all observed stars by 
counting them in strips (bins) parallel to the MS, as described 
in § 16.11 As for those distributions, the color index per bin is 
the average, corresponding to that at magnitude TO814 — 23. 
The distributions of Fig. [19] show an extraordinary similarity 
for the cases of LH 60, LH 63 and LH 72. Indeed, the fit- 
ting coefficients of the best-fitted Gaussian functions (drawn 
in red) are almost identical. Specifically, the peaks /i of the 
fitting functions are equal to 1.32, 1.29 and 1.41, and the cor- 
responding widths cr are 0.23, 0.23 and 0.25 forLH 60, LH 63 
and LH 72 respectively. The case of LH 88 is again an outlier 
and it is shown only for completeness. The CMD distribu- 
tion of the 'best candidate' stellar members for LH 88 is not 
represented by a single but by two Gaussians, with the first 
component of the fit (drawn in blue) being very narrow. This 
clearly suggests that the stars flagged in the previous section 
as best candidates for members of the star-forming region, and 
especially those at bluer colors, are in fact misidentified red- 
dened MS stars. The smaller red Gaussian component may 
represent traces of the PMS population of the system, which 
is contaminated by reddened evolved stars, but this result can- 
not be conclusive. 

The CMD distributions of the stellar members of LH 60, 
LH 63 and LH 72 show that these PMS stars cover a broad 
area of the CMDs of these regions. Indeed, the loci of PMS 
stars in the CMDs of Milky Way star-forming regions often 
show a wid ening, which could be evidence for age-spread 
jBriceiio eT al. 2007). This is also observed in young clus- 
ters and associations of both the Magellani c Clouds imaged 
with HST (see e.g., iGouliermis et al.ir2010l) . and our results 
so far show that the regions investigated here are not ex- 
ceptions. However, simulations performed previously by us 
have shown that characteristics of PMS stars, such as variabil- 
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TABLE 6 

Coefficients of the best-fitting Gaussian component 2 of the fits to each of the individual 
cross-sections distributions of faint stars across the cmds. 
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FWHM2 
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FWHM2 




M2 


FWHM2 




/^2 


FWHM2 
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17 


0.79 


0.49 


11 


0.85 


0.45 
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1.00 


0.58 
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1.17 


0.43 
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21 


0.93 


0.51 


14 


0.90 


0.45 


14 


0.92 


0.66 


14 


0.93 


0.45 
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53 


0.86 


0.41 


31 


0.89 


0.42 


14 


1.09 


0.76 


24 


0.98 


0.49 


4 


37 


1.09 


0.60 


24 


1.07 


0.70 


18 


1.17 


0.86 


18 


1.14 


0.53 


5 


31 


1.25 


0.76 


32 


1.13 


0.73 


20 


1.32 


1.05 


20 


1.09 


0.81 


6 


14 


1.45 


0.73 


16 


1.30 


0.89 


13 


1.46 


0.74 


6 


1.61 


0.14 



ity, binarity and circumstellar extinction can cause consider- 
able scatter of the positions of the PMS stars in the CMD, 
giving fals e evidence of a n age-spread within the systems 
jHennekem per et all l2008t iDa Rio. Gouliermis. & Gennarol 
1^10). Consequently we cannot attribute the observed CMD 
broadening of PMS stars solely to age-spread. 

Nevertheless, the distributions of Fig. [19] convey valuable 
information about the most probable ages of the systems and 
their indicative age-spreads. Taking the peaks of the cumula- 
tive distributions at 'face value' one can conclude that LH 63 
should be somewhat older than LH 60 and LH 72, because it 
corresponds to a bluer color, and thus to an older age. How- 
ever, considering the PMS stellar distributions in the individ- 
ual cross-sections across the CMDs for all three regions, they 
are all consistent with each other, peaking at colors along 
isochrones of ages between ^ 3 and 5 Myr, with the cases 
of LH 60 and LH 72 showing to be somewhat closer to the 
^ 3 Myr isochrone. As far as the widths of the PMS distribu- 
tions are concerned, they overlap with each other so well so 
that one cannot conclude any specific age difference among 
the three clusters. The derived age span of ^ 3 to 5 Myr is 
in excellent agreement with the age of the star-forming region 
LH 95, located at the north-western edge of the super-giant 
shell LMC 4, derived by a s elf-consistent age determination 
techn ique developed by us (IDa Rio. Gouliermis. & Gennarol 
120101) . In a subsequent study we will apply this technique to 
the data presented here, in order to constrain further the ages 
of the systems. It is worth noting that taking also the widths 
of the distributions at 'face value', they correspond to a sig- 
nificant spread between 2 and 10 Myr This result, however, 
is to be tested in our forthcoming study. 

Another factor to be considered as responsible for the CMD 
broadening of PMS stars is differential reddening. How- 
ever, while the PMS stars are expected to be somewhat dis- 
located due to reddening, our measurements of visual ex- 
tinction, discussed in § |3] show that differential reddening 
cannot either be fully responsible for the scattering of PMS 
stars in the CMDs of these three regions. Our measure- 
ments show that Ay in all investigated systems, except of 
LH 88, is quite low with a small variation with a maximum of 
Av — 0.68 mag (for LH 63), corresponding roughly to color 
excess of £;(F555W - F814W) ~ 0.27 mag, much smaller 
than the FWHM of the observed PMS broadening. On the 
other hand, LH 88 suffers from the highest extinction, with a 
maximum of Ay ~ 1.5 mag, which appears to be adequate to 
account for the observed CMD spread of stars. 

In conclusion, for the regions LH 60, LH 63 and 
LH 72, the cross-sections distributions of the PMS stars with 
771814 ^ 21.5 seem to be quite similar to each other in re- 
spect to their total numbers of stars, to their peaks and widths, 
as well as to the average color index at which the peak in stel- 
lar numbers appears. These extraordinary similarities clearly 



imply that the PMS stars found in three different young clus- 
ters, embedded in star-forming regions along the periphery of 
LMC 4, have possibly similar characteristics, and probably 
share a common star formation history. 

8. SUMMARY AND CONCLUDING REMARKS 

In this paper we present the first part of our investigation 
of PMS stellar populations in the LMC. Our targets of in- 
terest are four star-forming regions, LH 60, LH 63, LH 72 
and LH 88, located along the rim of the super-giant shell 
LMC 4. We present the reduction of the multi-wavelength im- 
ages taken with HST WFPC2 within our program GO- 1 1547 
(PI: D. Gouliermis), and our photometric analysis. We de- 
termine the accuracy and completeness of our photometry, 
and we measure the visual extinction towards the regions. 
We identify the stellar populations comprised in the observed 
fields in terms of the constructed CMDs. We evaluate the 
contribution by the field stellar populations with the use of 
archival WFPC2 images of the local LMC field. We make use 
of our rich photometric catalogs in the F555W and F814W 
(V- and /-equivalent) filters to assess the stellar content of 
the clusters embedded in the observed regions. 

We demonstrate that the true stellar content of the star- 
forming regions comprises both bright young stars observed 
at the UMS of their CMDs and faint stars still in the PMS 
phase of their evolution. An outlier from this general behav- 
ior is the region of LH 88, the high extinction of which does 
not allow a clear identification of its true populations. The 
identification of the PMS stars in the observed clusters is per- 
formed through the assessment of the contamination of the 
observed populations by the stellar content of the LMC field. 
This is achieved qualitatively by the constructio n of t he differ- 
ential Hess diagrams of the observed regions (ii l5.1b . and the 
derivation of the "clean" CMDs of the clusters after the statis- 
tical subtraction of the contaminant field population from the 
observed CMDs (^|5^. 

The probabilistic identification of the PMS stars in the ob- 
served regions is further performed, in a quantitative manner, 
by the distributions of the stars along cross-sections of the ob- 
served CMDs (§|6]l. We particularly focus on the faint PMS 
populations of the clusters, and we constrain our analysis to 
the fainter part of the CMDs with nisu ^ 21.5, which is less 
contaminated by evolved field populations. We construct the 
number distributions of the stars across each observed CMD 
along six selected magnitude ranges and we fit these distri- 
butions with a two Gaussian components. The first compo- 
nent represents the MS field stars observed in each region, 
while the second corresponds to the PMS member stars of 
the region. For each star we assign a membership probabil- 
ity, p, derived from the ratio of the second fitted Gaussian 
to the sum of the two components along each cross-section. 
We isolate the stars with p > 95%, as the best candidates of 
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being members of the star-forming regions. All these can- 
didates are PMS stars. The CMD distributions of these stars 
show an extraordinary similarity to one another for the regions 
LH 60, LH 63 and LH 72, suggesting similar characteristics. 
Considering that the peaks of these distributions represent the 
most probable ages of the regions, their similarity also sug- 
gests that all three regions may share a common recent star 
formation history. This result is quite important, since these 
regions are located at different parts of the boundaries of the 
super-giant shell LMC 4, and therefore our findings suggest 
that star formation around the shell may have occurred almost 
simultaneously. Nevertheless, these distributions show a def- 
inite widening of the loci of faint PMS stars in the observed 
y/-equivalent CMDs, which should be further investigated. 

The spread of PMS stars along the color axis in the CMDs 
of star-forming regions in the Magellanic Clouds is a well- 
documented phenomenon, which demonstrates that the par- 
ticular location of a PMS star in the CMD i s a very complex 
function of its intrinsic properties (see, e.g.. lGouliermis et alj 
|2010>) . Low-mass PMS stars in such regions, being the coun- 
terpart of Galactic T Tauri stars, suffer from rotational vari- 
ability, non-periodic variability due to accretion, circumstel- 
lar extinction, and binarity. These characteristics dislocate 
the stars from their original positions in the CMD in an 
unpredictable manner, which may be misinterpreted as an 
age-spread within the host stellar system. The only accu- 
rate approach to quantify this effect and to confirm or re- 
fute the spread in ages is through detailed simulations of 



synthetic CMDs and their comparison to the observed ones. 
We are currently undertaking such an investigation with the 
accurate statistical determination of the masses and ages of 
the identified PMS stars through t wo self-consistent tech- 
niques established by us (see Da Rio. Goulier mis. & Henninp 
2009; Da Rio, Gouliermis, & Gennaro 2010). Our subse- 
quent investigation of the observed star-forming regions will 
include the multi-wavelength characteriza tion of the bright 
stella r content of the clusters (see, e.g., iRomaniello et alj 
l2002h . the determination of circumstellar accretion of bright 
PMS sta rs through their U- and Ha-excess emission (e.g., 
iDe Marciii. Panagia. & Romanienoll2010 ). and the identifica- 
tion of embedded massive candidate young stellar objects and 
ultra-compact Hll regions via the syner gy of our observa tions 
with those in longer wavelengths (e.g., IGruendl & Chul2009i) . 
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